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Abstract 

Single pion production off nucleons is studied in the framework of relativistic baryon 
chiral perturbation theory at tree level with the inclusion of the terms from the 
dimension two effective pion-nucleon Lagrangian. The five appearing low-energy 
constants are fixed from pion-nucleon scattering data. Despite the simplicity of 
the approach, most of the existing data for total and differential cross sections as 
well as for the angular correlation functions for incoming pion kinetic energies up 
to 400 MeV can be satisfactorily described. 



PACS nos.: 25.80.Hp , 12.39.Fe , ll.30.Rd 



1 Introduction 



Single pion production off nucleons has been at the center of numerous experimental 
and theoretical investigations since many years. One of the original motivations of these 
works was the observation that the elusive pion-pion threshold S-wave interaction can be 
deduced from the pion-pole graph contribution. A whole series of precision experiments 
at PSI, TRIUMF and CERN has been performed over the last decade and there is still 
on-going activity. On the theoretical side, chiral perturbation theory has emerged as a 
precision tool in low energy hadron physics. Beringer considered the reaction vrA^ ttttN 
to lowest order in chiral perturbation theory []T|. Low-energy theorems for the threshold 
amplitudes Di and -02^3 w^re derived in 0. These are free of unknown parameters 
and not sensitive to the tttt interaction beyond tree level. A direct comparison with 
the threshold data for the channel n^p 7r~^TT^n, which is only sensitive to Di, lead 
to a very satisfactory description whereas in case of the process 7T~p vr^vr'^n, which 
is only sensitive to D2, sizeable deviations were found for the total cross sections near 
threshold. These were originally attributed to the strong pionic final-state interactions in 
the I-jtiT = channel. However, this conjecture turned out to be incorrect when a complete 
higher order calculation of the threshold amplitudes Di^2 was performed 0] . In that paper, 
we investigated the relation between the threshold amplitudes Di and D2 for the reaction 
TcN —>■ ttttN and the tttt S-wave scattering lengths Oq and ag in the framework of heavy 
baryon chiral perturbation theory to second order in small momenta, i.e. the pion mass 
(which is the only small parameter at threshold). The pion loop and pionic counterterm 
corrections only start contributing to the ttttN threshold amplitudes at second order. One 
of these counterterms, proportional to the low-energy constant £3, eventually allows to 
measure the scalar quark condensate, i.e. the strength of the spontaneous chiral symmetry 
breaking in QCD. However, at that order, the largest contribution stems indeed from 
insertions of the dimension two chiral pion-nucleon Lagrangian, which is characterized by 
a few coupling constants called Cj. In particular this is the case for the amplitude D2. To 
be specific, consider the threshold amplitudes D12 calculated from the relativistic Born 
graphs (with lowest order vertices) and the relativistic Q-terms expanded to second order 
in the pion mass. This gives 

^Born^^c, ^ (2.33 + 0.24 ± 0. 10) fm^ = (2.57 ± 0. 10) fm^ , (1) 
^Born ^ jjc, ^ ^-Q,Qi _ 2.85 ± 0.06) fm^ = (-9.46 ± 0.06) fm^ , (2) 

which are within 15% and 5% off the empirical values, Df^^ = (2.26 ± 0.10) fm^ and 
D^^^ = (—9.05 ± 0.36) fm^, respectively. It appears therefore natural to extend the same 
calculation above threshold and to compare to the large body of data for the various 
reaction channels that exist. It was already shown by Beringer [|l| that taking simply the 

#^These are related to the more commonly used Aw and A32 by -432 = VlODi and = — 21)1 — 3I?2- 
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relativistic Bom terms does indeed not suffice to describe the total cross section data for 
incoming pion energies up to 400 MeV in most channels. Such a failure can be expected 
from the threshold expansion of D2-, where the Born terms only amount to 75% of the 
empirical value. We therefore expect that the inclusion of the dimension two operators, 
which clearly improves the prediction for D2 at threshold, will lead to a better description 
of the above threshold data. In particular, it will tell to which extent loop effects are 
necessary (and thus testing the sensitivity to the pion-pion interaction beyond tree level) 
and to which extent one has to incorporate resonance degrees of freedom like the Roper 
and the A-isobar as well as heavier mesons (ex, p, uj) as dynamical degrees of freedom (as 
it is done in many models, see e.g. [Q). Our calculation is in spirit closest to the one 
of ref.|^, in which Beringer's Born terms where supplemented by explicit A and Roper 
(tree) contributions. Clearly, the inclusion of the resonances as done in that paper is not 
based on a consistent power counting scheme but rather it is argued that phenomenology 
demands the extension of the effective Lagrangian to include these higher mass states. It 
will be of particular importance to directly compare our results to the ones presented there 
for the abovementioned reasons. Since we will only calculate tree level diagrams with at 
most one insertion from the dimension two pion-nucleon Lagrangian, it is advantageous to 
treat the nucleons as fully relativistic Dirac fields, since then one automatically generates 
all 1/m suppressed terms with fixed coefficients. Here, m denotes the nucleon mass. Only 
at next order, when loop graphs have to be included, the problems due to the additional 
large mass scale given by the nucleon mass appear . Therefore, the calculation presented 
here can not easily be extended to higher orders in the chiral expansion. At this point, 
we stress that the amplitudes calculated from the tree graphs as done here are, of course, 
purely real. How severe this approximation is can only be judged when a full one loop 
calculation above threshold has been performed. 

This paper is organized as follows. Section 2 contains some formal aspects including the 
definition of the T-matrix and of the pertinent observables. In section 3 we briefly review 
the effective Lagrangian underlying the calculation and discuss the tree level diagrams 
to be evaluated. Section 4 contains the results and discussions thereof. The appendices 
include more details about the kinematics, some lengthy formula appearing in the expres- 
sions for the observables and the explicit expressions for the invariant amplitudes. 

2 Single pion production: Formal aspects 

In this section, we outline the basic framework concerning the reaction vrA^ nirN above 
threshold. Some of the lengthy formulae are relegated to the appendices. 
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2.1 T-matrix and invariant amplitudes 

We seek the T-matrix element for the process 

7r«(fc) + N{p^) ^ n\q,) + n^q^) + N{p^) , (3) 

with denoting a nucleon (neutron or proton) and tt" a pion of (cartesian) isospin a. 
This process is characterized by five independent Mandelstam variables, 

s = (p, + A;)2, si = (gi+p2)', S2 = {q2+P2T, h = {qi-kf, t2 = {q2 - kf . {4) 

All scalar products between the various momenta can be expressed in terms of these, see 
appendix A, and at threshold (gi — q2 — 0, in the center-of-mass frame), we have 

= (m + 2M.)^ sr = sr = {m+M^r, tr=tr^-Mf^^^^^±^. (5) 

m + zM^ 

The nN — > nnN transition matrix element can be expressed in terms of four invariant 
amplitudes, denoted Fj(s, si, S2, ^i, ^2) — 1, 2, 3, 4), 

T = i U2i^[F^ + (?i + (^2)F2 + (gl - (^2)^3 + {(^l(^2 - (/2(/i)F4ui . (6) 

with ui^2 the Dirac spinor for the in- and outgoing nucleon, respectively. In the complete 
relativistic tree calculation to order g^, these are rational real functions of the five Man- 
delstam variables. The isospin decomposition of the invariant amplitudes reads (pulling 
out a common prefactor composed of coupling constants) 



9A 

8F3 



(7) 



with gA the nucleon axial-vector coupling constant and F^^ the pion decay constant, 
respectively. The crossing properties of the isospin amplitudes can be readily deduced. 
Crossing the two out-going pion lines, {qi,b) <-> (^2, c), amounts to 

Bj{s,Si,S2,ti,t2) ^ eiejBl{s,S2,Si,t2,ti) , ei,2 = +1 , £3,4 = -1 ■ (8) 

In the physical reaction channels the relevant amplitudes are given by 

F^-§-,tvA, (9) 

j = l 

with r)j channel dependent isospin factors. We calculate the amplitudes Bj in the isospin 
limit with the charged pion mass (M^r = 139.57 MeV) and the proton mass, m — 
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938.27 MeV, in order. Isospin breaking is done by shifting the kinetic energy of the 
incoming pion, T^, from the isospin symmetric threshold 



3M^ 
2m 



170.71 MeV 



(10) 



to the correct threshold. The pertinent isospin coefficients rjj and numerical values for 
this shift, called ST^^, for the five experimentally accessible reaction channels are 

, = 2a/2 , 7^3 = , 7^4 = , = +1.68 MeV (11) 

, 7^2 = 1 , 773 = -1 , 7/4 = 1 , = -5.95 MeV (12) 

V2,7]2 = V2,r]3 = 72,7/4 = 0,5T, = +1.68 MeV (13) 

V2,7]2 = 0,773 = 0,774 = 0,6T^ = -10.21MeV (14) 

, 7/2 = 1 , 7/3 = 1 , 7/4 = 1 , = -5.95 MeV . (15) 



tt'^P — 


> Tr'^TT'^n : 


^1 


tt'^P - 


-> n'^n^p : 


^1 


Tr~p — 


> Tr'^Tt~n : 


^1 


7l~p - 


-^nW : 


^1 


n^p - 


n^n^p : 


^1 



2.2 Observables 

The total cross section is a four-dimensional integral over a quadratic form in the ampli- 
tudes Fj, 



C^tot(^7r) 



(47r)47n^7;(7; + 2M,) 



/ / dujiduj2 



dxi 

1 JQ 



. ^ y,,F:F, . (16) 



The weight functions yij are given in appendix A, and 5 is a Bose symmetry factor, 
iS = 1/2 for identical pions in the final state and 5 = 1 otherwise. The center-of-mass 
(CMS) kinematics is 







s-m^ + Ml 



\k\ 



ki-Ml, 



Si ^ s - 2^/s uJs-i + , U = 2{MI - koUi + \k\\qi\xi) {i 

X2 — XiZ + ^(1 — x\){l — Z^) COS (j) , 

z l?i||?2| = ^1^2 — \/s{<^\ -\-0J2) + Ml + ^(s - m^) . 



1,2) , 



(17) 
(18) 

(19) 
(20) 



Here, is the (auxiliary) angle between the planes spanned by q\ and k as well as q\ and 
q2. In accordance with the experimentalists convention, we have chosen the coordinate 
frame such that the incoming pion momentum k defines the z-direction, whereas gl, the 
momentum of tt*" lies in the xz-plane. The polar angles Q\ and Q2 of the outgoing pions 
(with Xi — cos Qi) are in general non-vanishing and so is the azimuthal angle 02 of tt'^- By 
construction, the azimuthal angle of 'k^ is zero. 



5 



The double differential cross section in the CMS depends on T^, tui, and 6i. For this more 
exclusive quantity the energy ui and the polar angle 9i of vr* are detected, 



2„ OC p-K 4 

dU2 / d<P J2 yV^:Fj (21) 



with the boundaries of the three-body phase space given by 

1 



^2 



2{s - 2^suJi + Ml) 



(v^ - uji){s - 2^uji -m^ + 2Ml) 



[22) 



Also measured are triple differential cross sections. These depend on Tj,, uji, 6i, 62, and 
02- The two angles of tt*^ are measured in addition, thus the three-particle final state is 
kinematically completely determined, 

d^a \qi\\q2\S A ~ , . 

1^ VijFi Fj , (23) 



dujidnidn2 {A-n)^^\k\E2 CJ^i 
with s, si, S2, ti, t2 as before and 

f, ^ ^^A ^2(1(3 - m^) - y/j^l) + M^juJi +UJ2-V^) 

= + = ^F^M^ ' ^'^^ 

with the final state nucleon energy. The cosine of the angle between qi and q2 is given 

by 

z = cos 61 COS 62 + sin 61 sin 6'2 cos </)2 , (25) 
and it is needed to express the energy of vr"^ as 

1 



002 



2[{^s - iu^)^ - z^\q,\ 



-|(v^-cji)(s-2v^cui -m2 + 2M^) 



-z\qi 



^{s- 2^/^uJl - m2)2 - 4M2(m2 + (1 - z^)\qi\'^)^ . (26) 

From the triple and double differential cross sections, one defines the angular correlation 
function W at fixed beam energy T^, 

WM,O.M = *-^;j^J^^. (27) 
This completes the necessary formalism for our calculation. 



3 Calculation of the tree level amplitudes 

In this section, we briefly discuss the effective chiral Lagrangian underlying our calcula- 
tions. For more details we refer to the review and to the paper 0. We then show and 
discuss the various classes of tree diagrams contributing to the process vrA^ — > nnN to 
second order in the chiral expansion. 
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3.1 Effective Lagrangian 

At low energies, the relevant degrees of freedom are hadrons, in particular the Goldstone 
bosons linked to the spontaneous chiral symmetry breaking. We consider here the two 
flavor case and thus deal with the triplet of pions, collected in the matrix U{x) = u'^{x). 
It is straightforward to build an effective Lagrangian to describe their interactions, called 
^TT^. This Lagrangian admits a dual expansion in small (external) momenta and quark 
(meson) masses as detailed below. Matter fields such as nucleons can also be included 
in the effective field theory based on the familiar notions of non-linearly realized chiral 
symmetry. The pertinent effective Lagrangian is called C-j^n, consisting of terms with 
exactly one nucleon in the initial and the final state. The various terms contributing to 
a process under consideration are organized according to their chiral dimension, which 
counts the number of derivatives and/or meson mass insertions. Here, we work to second 
order in the corresponding small parameter q (which is a generic symbol for an external 
momentum or pion mass). Consequently, the effective Lagrangian consists of the following 
pieces: 

£eff = c^:! + /:a + /:a , (28) 

where the index (z) gives the chiral dimension. The form of C^^J + C^^lf is standard. Let us 

(2) 

discuss in some more detail the terms appearing in C;^. For the reasons outlined before, 
we treat the spin-1/2 fields, i.e. the nucleons, relativistically. The terms of the dimension 
two effective relativistic pion-nucleon Lagrangian relevant to our calculation read 0] 

£S = Ci^'^Tr(x+) + ^^tlf. Du niu'^u") - B^B. ^Triu^^u") 

+ Cs^uX^ + c^-^a^,{u^, , (29) 
where \l/ denotes the relativistic nucleon field and with 

X+ = X"^^ + ) = 2 f/ . (30) 

The quantity x contains the light quark mass m and external scalar and pseudoscalar 
sources |]^ and is the standard chiral covariant derivative acting on the pion fields. 
The low-energy constants Cj are normalized such that we can identifiy them with the 
corresponding low-energy constants of the heavy baryon Lagrangian (truncated at order 
g^) (for definitions, see e.g. 0). Their numerical values have been determined in from 
a tree level fit to pion-nucleon scattering (sub)threshold parameters, 

ci = -0.64 ±0.14, 4 = -5.63 ±0.10 , 4' = 7.41 ± 0.10 , 

cg = -3.90 ±0.09, C4 = 2.25 ±0.09 , (31) 

with all numbers given in GeV^^. It is important to stress that the effective Lagrangian 
we use does not incorporate any resonances like e.g. the A(1232) or the iV*(1440) as 
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dynamical degrees of freedom. Their contribution is encoded in the numerical values of 
the low-energy constants q, for a more detailed discussion see [|T^. Clearly, one might 
doubt that above threshold such a procedure is sufficient, since for example it does not 
account for the effects related to the widths of these resonances. If one expands the 
resonance propagators in powers of momentum transfer over resonance mass, one gets to 
leading order such dimension two pion-nucleon couplings as used here (plus an infinite 
series of higher order terms). Our study should therefore reveal at which kinematics and 
for which observables explicit resonance effects can be most easily seen. This would be 
helpful for pinning down various resonance couplings which are not too well constraint by 
the present models (when global fits are performed). Finally, we stress that the dimension 
two operators encode also information about t-channel meson exchanges, e.g. the constant 
Ci is essentially saturated by scalar meson (correlated 27r, to be precise) exchange and C4 
receives an important contribution from the p UTOf . 

3.2 Chiral expansion of the invariant amplitudes 

Let B he a. generic symbol for any one of the Bl. The chiral expansion considered here 
takes the form 

Consider first the Born terms with insertions solely from d]^}^ shown in Fig. 1. Diagram A 
involves the A^A^Svr-vertex genuine to ttN —>■ tc-kN at tree level. The lowest order four- 
pion interaction appears in the pion pole graph (B). The graphs of the classes C, D, E and 
F contain one pion-nucleon vertex and one Weinberg NNtttt interaction, thus these scale 
as qa- The graphs belonging to the classes G, H and I each contain three pion-nucleon 
vertices and thus scale as g\. In Fig. 2, we show the term with exactly one insertion from 

(2) 

£^jv, i.e. these are proportional to the low-energy constants q. Note that there is no 
higher order A^A^Svr-vertex from the dimension two Lagrangian because of isospin. Note 
also that the graphs of type J, K, L and M are topologically identical to the ones of classes 
N, O, P and Q. However, while the operator proportional to C4 contributes to the Bf, the 
ones ~ ci,C2,C2 and C3 can not. The explicit formula for the Bf for the various (classes 
of) diagrams are given in appendix B. 

4 Results and discussion 

In this section, we show the results for the total and double differential cross sections as 
well as the angular correlation function in comparison to the available data. Our input 
parameters are = 139.57MeV, m = 938.27MeV, = 92.4 MeV and gA = 1.32 as 
determined from the Goldberger-Treiman relation qa = Ft^Qt^n/^ with Qt^n = 13.4. For 
the low-energy constants Cj, we take the central values given in Eq. (|3lD . Setting Cj = 
corresponds to the tree level (dimension one) calculation of Beringer |^ . 
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4.1 Total cross sections 



First, we discuss the data from threshold up to incoming pion kinetic energy of = 
400 MeV and then take a closer look at the threshold regions, T,r < 210 MeV. In all 
graphs, the solid line gives the result with the q as given in Eq.(^) and the dashed line 
the one with the q = 0. 

TT^p TT+TT+n, Fig. 3: The data are from refs. |ll| ||12| [|1^ The lowest order cal- 



culation does already quite well. Inclusion of the dimension two terms leads to a too 
large cross section. This can be traced back to the fact that the Born term result for 
the threshold amplitude Di is closer to the empirical number than the one with the Cj 
corrections, compare Eq.(|l]). 

TT+p — > TT+TT^p, Fig. 4: There are recent data from [TH] and a few older data in this kine- 



matical range, see refs. ||T6| While the lowest order calculation tends to under- 

estimate the data above the threshold region, the q corrections lead to a somewhat too 



large total cross section, in particular when compared only to the data of |T^ . 



TT p 



Fig. 5: The data are from refs. HI O [P H H H [H H 13- The 



lowest order calculation clearly underestimate the cross section even close to threshold. 
Inclusion of the dimension two terms leads to an astonishingly good description of the 
total cross section up to T^r = 400 MeV. 



IT p 



Fig. 6: The data are from refs. [^] The lowest order 



calculation clearly underestimates the cross section even close to threshold. Inclusion 
of the dimension two terms leads to a good description of the total cross section up to 
= 230 MeV and a slight underestimation for larger pion energies. 



71^ p 71^71^ p, Fig. 7: The data are from refs. ||2T| |^. Inclusion of the dimension 
two terms leads to an improved description of the total cross section up to = 250 MeV 
and a slight overestimation for larger pion energies whereas the lowest order result is 
consistently below the data for > 230 MeV. 

Threshold region: The threshold regions for all five channels are shown in Figs. 8 and 9. 
In all but one case, i.e. for 7i~^p —>■ 7c^7c~^n, the inclusion of the terms ~ q clearly improves 
the description of the data. As already noted in the introduction, this improvement is 
not due to the strong pionic rescattering in the isospin zero, S-wave as originally thought 
0. It should be noted, however, that only for 7c~^p Ti^Ti^n and n^p n^Tc^n, precise 
data close to threshold exist 



4.2 Double differential cross sections 

Double differential cross sections have been published by Manley [^] and by Miiller et al. 



3^ for the reaction tx p ^ tt+tt n. In Fig. 10, we show the results in comparison to the 

1.242 GeV and for = ui - = 6.5, 10.95 and 15.35 MeV, in 

1.262 GeV and = cui - = 10.4, 17.5 



data of at 

order. Fig. 11 comprises the results for ^/s 
and 24.6 MeV. In all cases, our full calculation describes well the data whereas the leading 
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order calculation is systematically below these. However, we are not able to describe the 
data of 0] at y/s = 1.301 GeV as shown in Fig. 12. The two solid lines in that figure refer 



to the angular range of Xi = cos 6*1, 0.03 < Xi < 0.35. The strong peaking of the data 
around Tt^+ = uji — M^^ ~ 25 MeV can not be reproduced. Two comments are in order. 



The experiment of |^ was not intended primarily to measure the double differential cross 
section but rather the angular correlation function, see the next paragraph. It is, however, 
claimed in |Q that their results for d?a / duoidVti agree with the ones of Manley |Q |3^ 



once a scaling factor to account for the different incident pion kinetic energies is applied 
and the different laboratory angles of both experiments are considered. Other double 
differential cross sections for the same process measured by the Erlangen group are given 
in 1^ These span the range of the incoming pion kinetic energy T^- from 218 to 330 
MeV. The data for the lower energies are well described within our approach, as shown in 
Fig. 13. For beam energies T^- = 246 and 330 MeV, our results underestimate the data, 
but not as strongly as for the data of |3^ at T^^- = 284 MeV. 



4.3 Angular correlation functions 

An extensive study of the angular correlation function W for tt^p -K^-K^n at ^/s 
1.301 GeV is given in [^. For six values of the momentum of the positive pion, \qi \ 



60.5, 86, 100, 113, 127 and 151.5 MeVp^ tables are given for W a.s a. function of 02 (with 



6*2 fixed) and of 62 (with 02 fixed). We scanned through all of these data and found the 
following results. For fixed 6'2, the angular correlation function is well described as long as 
02 > 50°. For smaller azimuthal angles, we consistently underestimate the data, a typical 
example is shown in Fig. 14. On the other hand, for the fixed values of 02 given in |Q , 



the 62 dependence is satisfactorily described, two examples are given in Fig. 15. We have 
also considered the angular correlation functions given in and ||3^ and found similar 



trends, i.e. in general a good description of W as long as the dimension two terms are 
taken into account. 

4.4 Comparison to other calculations 

As already stated, the calculation of is in spirit closest to ours. They use the relativistic 
pion-nucleon Lagrangian to leading ordeiP^ and add explicit A and Roper degrees of 



freedom, but no heavy meson (a, p) exchanges. Overall, their results are not significantly 
better than ours. It is important to note that these authors did not try a best fit by fine 
tuning the various coupling constants. Our conclusion is similar to theirs, namely that for 
really testing the constraints from chiral symmetry, one needs more precise data at low 
energies. The model of @] is more elaborate concerning the inclusion of meson exchanges 



"^^We use the central values of the pertinent momentum bins. 
#3The are not using heavy baryon CHPT as claimed in their paper. 
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and is more tuned to the resonance than to the threshold region. Similarly, the work 
of the Erlangen group ^ shows the importance of the non-linear pion couplings at tree 
level but does not include higher derivative operators as done here. Our appraoch lends 
more credit to the conclusions of these papers that the intricate dynamics of the reaction 
vrA^ —>■ tcttN reveals itself best in double differential cross sections and angular correlation 
functions but not in the total cross sections. 



5 Summary and conclusions 

In this paper, we have studied the reaction vrA^ ttttN in relativistic baryon chiral 
perturbation theory at tree level including the dimension two pion-nucleon Lagrangian. 
This study was motivated by the fact that at threshold the dominant corrections to the 
lowest order tree graphs indeed come from tree graphs with exactly one insertion from £^j(f 
and not from loop graphs, as shown in . The values of the relevant low-energy constants 
can be fixed from a few pion-nucleon scattering data. These low-energy constants encode 
information about baryonic and mesonic resonance excitations, in particular related to the 
A, the Roper, the p and the correlated two-pion exchange in the scalar-isoscalar channel 
|TD|. Consequently, all our results concerning single pion production are based on a truely 



parameter-free calculation. The pertinent results of this study can be summarized as 
follows: 

o For pion energies up to T^^ = 250 MeV, in all but one case the inclusion of the 
contribution ~ Cj clearly improves the description of the total cross sections, most 
notably in the threshold region for tt^p —>■ TT^n'^n. Up to T^r = 400 , the trend of the 
data can be described although some discrepancies particularly towards the higher 
energies persist, see Figs. 3. . .9. This indicates that our approach is best suited for 
the threshold regions. 



o 



Double differential cross sections for 7r~p Ti^Ti~n at incident pion energies below 
= 250 MeV are well described (see Figs. 10,11,13). We are not able to reproduce 
the strong peak of d?a / duidVLi as a function the outgoing tt"*" energy at y/s = 
1.301 GeV found in [H, compare Fig. 12. 

o The angular correlation function for the reaction 7T~p Tr+vr^n as given in ||3^ can 
be reproduced fairly well for all energy bins and sets of polar and azimuthal angles 
given in that paper, with the exception of small polar angles for the final-state 
negative pion, compare Figs. 14 and 15. 
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A Kinematics and weight functions 

First, we express all possible scalar products of the nucleon and pion momenta in the 
initial and final state in terms of the five Mandelstam invariants, see Eq.(D, and the 
particle masses, 

2p^.k = s-m^ - Ml, 2pi-qi = s-S2 + h- M^, 2pi ■ q2 = s - si +t2 - M^, (A.l) 

2pi-p2 = Si + S2-s-ti-t2 + m^ + Ml, 2k-qi = 2Ml-ti, 2fc ■ = 21/^ -ta, (A.2) 

2k-p2 = s + ti + t2-m^ - 3Ml, 2qi ■ q2 = s - Si - S2 + (A.3) 

2gi ■P2 = si-m^ - Ml, 2q2 ■ P2 = S2 - - Ml . (A.4) 
The weight functions yij appearing in Eqs.(|l^,^,^) are given by: 

yii = si + S2 - h - t2 - s - + Ml , 
yi2 = m{2{si + S2 - s - m^) - ti - 12) , 
yi3 = rn{t2 - ti) , 

Via = {si- S2){s - si- S2 + m^) + (ti-t2){m^ + Ml) + sit2- S2ti, 
y22 = isi + S2){^m^-Ml) + {ti+t2)is-m^) + s^-3s{2m^ + Ml) 



-3m^(m'^ + Ml) + 2M^ 



1/23 = {si - S2) {s-m^- 2Ml) + (t2 -ti){m^ + Ml) + Sih - ^2^2 , 

?/24 = m{2si-2s2 + ti-t2){s - Si- S2 + m^ + 2Ml) , 

1/33 = {si + S2){2s-2m^-3Ml) + {ti+t2){2Ml-m^-s)-4siS2 

+2siti + 2s2t2 -s^ + s(2m2 + 3M^) + 3m^ - Sm^M^ - 2M^ , 
yu = m{2s + ti+t2-2m^ -4Ml){s - si- S2 + m^ -2Ml) , 
1/44 = {s - Si - S2 + rri^)[{s - m^){s + ti +t2 - m^) + {si - S2){s2 - Si +t2 - ti)] 

+Ml{4Ml{2s - Ml) + m^[3m^ - 8Ml + Q{s - si - S2) + 2{ti + ta)] 

-5S^ + 2s{Si + S2-ti- t2) -sj-sl + 6S1S2 + 2{Si - S2)it2 - ti)} . (A.5) 

and they are symmetric under interchange of the indices, yij = yji. 

B Invariant amplitudes Bf 

In this appendix we give the invariant amplitudes Bj for individual (classes of) diagrams. 
The Bj not specified are zero for the class of diagrams just considered. The following 
abbreviations for kinematical quantities are used: 

di = si — , 8.2 = S2 — 1 = 2Ml — s + Si — t2 , (B-l) 
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di = 2Ml - s + S2 - ti , d5 = m^ + AMI - ^ - - ^2 , 4 = s - . (B.2) 
Diagram A: 

^2 = 2, Bl^Am, Bl^Bl^-l. (B.3) 



Diagram B: 



1 _ Sm{di + d2-de- Ml) , ^ 4m(2M^ - ti - ta) „3 ^ 4m(ti - ts) 

-"1 O 71 J J J ) -"1 o 71 ^9 J J J ' 1 



2M2-di-d2-4 2M2-di-d2-4 2M2-di-4-4 

(B.4) 

Here, the separation of the invariant amphtudes between diagram A and B results from 
using the cr-model gauge, U — ^1 — t^/F^ + ir ■ tt/Ftt. We have checked that other 
parametrizations lead to the same result for diagram A plus B. 
Diagram C: 

El = Bl = 2m(s2 - si)dl^ , E| = 5^ = 1 , Bf = = -2mde ^ . (B.5) 
Diagram D: 

Bf = -Bl = 2m(4 - di)d^^ , 5| = -B^ = 1 , 5| = -5^ = 2md^^ . (B.6) 
Diagrams E: 

Bl = -2S? = 2m(rfi + d2- de) {d^^ + d^^) , Bl = -2S| = -3 , (B.7) 

Bl = -2Bl = 2m(ci2'^ - (^r^) , = 2S? = 2m{di + d2 - 4) ((^r^ - ^^2"^) , (B.8) 

= 2Sf = 1 , = 2S| = -2m(di^ + d^^) . (B.9) 

Diagrams F: 

Bl = -2^2 = 2m(cii + ci2 - 4)((i3 ^ + d^^) , = -2^^ = -3 , (B.IO) 

= -2S| = 2m(rfJ^ - c^3 ^) , Bf = -2Bf = 2m{di + d2 - de){d^^ - d^^) , (B.ll) 

= -25^ = -1 , 5^ = -25| = -2m{dl^ + ^) . (B.12) 

Diagrams G: 

Bl = g\mdl^{{s + 3m2)[(s - 81)^2 ^ + (s - S2)rfr^] - 4(s + m^)} , (B.13) 

Bl = = g\mdl\{s + 3m2)[(s - S2)cir' - {s - si)d:,^] + 2(si - S2)} , (B.14) 
^2 = gW^rn'^d^^ + 1 - 2m^(dr' + ^^2 ')] , Bl = B^ = = 2^im2(d2-i - dr') , (B.15) 

S| = Bt = gl[-l - 2m\dT' + d^')] , Bj = glm{s + 2,m^)d^\dl^ - ^2"^) , (B.16) 

Bl = Bl = g\mdl\2 + (s + ?>m^){d'{^ + ^2 ')] • (B.17) 
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Diagrams H: 

Bl = -Bf = g\m{2 + {s - S2)d^^ + {s - si)d^^ 

+ (de -di- d2)[(3m' + si)(ir'o?3"' + (3m' + S2)d^^d^^]} , (B.18) 

Bf = g\m{ {s — S2)di^ ~{^~ •si)c?2 ^ + {de — di — 0^2) [(3m' + Si)di^d^^ — (3m' + S2)d2^d^^} , 

(B.19) 

Bl = -Bl = 2gl[l + m'idT' + d^' + d^' + d^')] , B^ = 2^>'(dr' - d^' + d^' - d^) , 

(B.20) 

= -Bl = 2^>'(cir' - d^' - d^' + dt) , (B.21) 
B', = 2gWidT' + d^' - d^' - df) , (B.22) 
Bl = -Bl = g\m[{Am^ + di + d3)c^r^t^i"^ - (4m' + ^2 + ^4)^2^^^] , (B.23) 
54 ^ ^2^[(4^2 + + d^)d^^d^^ + (4m' + (^2 + (i4)c^2 ^C^] • (B.24) 

Diagrams I: 

B\ = ^^mds ^{(4m' + d5)[(4 - (^4)^3' + (4 - ^3)^1'] - 4(^5 + 2m')} , (B.25) 

Bl = -B{ = g\mdl\{d^ + 4m') [(4 - d^)dl^ - (rfg - 4)<3 ^] + 2(^4 - 4)} , (B.26) 
Bl =gl[l + 8m^d^' - 2m\d^' + d^^)] , Bl = -E'l = Bl = 2g\m\d^' - d^') , (B.27) 
= -Bl = g\[-l - 2m'(d3 1 + d^^)] , Bl = g\m{d, + Am')d^\d^^ - d^) , (B.28) 
Bl = -Bl = ^^(is ^ [-2 - (4 + 4m') (^3"^ + t^j')] ■ (B.29) 

Diagrams J: 

Bl = i{4(3m' + s) [AciMl + 03(4 - ^2 - cii)] + [^e - (^1 - ^2)'] (2^ + 4'(4/(2m') + 2)]} , 

(B.30) 

Bl = -4(4 + C2)/m - mdg ^[32ciM^ - 4(4 - 4)7^^ + 8^3(4 - 4 - 4)] , (B.31) 

^3 = 4(^1 - S2)/m , Bl = 24(s2 - si)c?6 ' ■ (B-32) 

Diagrams K: 

= 4-^{4(4m' + 4)[4ciM,' + C3(4-4-4)] + [rf5-(^^3-4)'](24 + 4(4/(2m') + 2)]}, 

(B.33) 

Bl = -4(4 + 4)/m - mdg ^[32ciM^ - 4(4 - 4)'/m' + 803(4 - 4 - 4)] , (B.34) 
Bl = 4 (4 - 4) /m , Bl = 24 (4 - 4)(i5 ^ . (B.35) 

Diagrams L: 

Bl = i6ciM^-2c3(4 + 4) + 4(4-4-4)[(4-4)4^ + (4-4)4^] 

+4/(2m')[(4 -Ml + ti/2)(4 - s + S2) + (4 - + ^2/2) (4 - s + si)] , 

(B.36) 
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Bl = 2cs{t2-h)+c',{de-di-d2){[2{M^-d,)-t2]di'-[2{M^-de)-h]d^'} 
+4/m'^[m\s2 - si + (t2 - ti)/2) + {s - Ml){d^ - d^) 

+t2/2(si - t2/2) - h/2{s2 - h/2)] , (B.37) 

Bi = 8c,M^m{d^' + d2^) + 2c3m[(t2 - 2M2)rff ^ + {t^ - 2M^)d^^] 

+c'2/m[m^ + - 3s + Si + S2 - 3/4(ti + ^2)] + c'27m[2(4 - M^) 
+(^1 + ^2)72 - (s + ^2/2 - - M2)2(i^i - (s + ti/2 -m^- Mlfd2^] , 

(B.38) 

Bl = SciMX^r' - (^2"') + 2c3m[(i2 - 2M2)(i-i - {t^ - 2M^)d^^] 
+4/m[si - S2 + 3/4(ii - ^2)] + 4'M[(^2 - h)/2 

-{s + t2/2 - - M^fd^' + {s + h/2 - - M^)^^-^] , (B.39) 

52 = Bl + 4(4 - c^s) /m , (B.40) 

Bl^Bl + 4 (3s - - si - S2 - 4) , (B.41) 

Bl = 4[(2(M2 - 4) - t2)rfr' - (2(M2 - d,) - h)d2^] , (B.42) 

Bl = 4[2 + (2(M,2 - 4) - ^2)^:' + (2(M,2 - 4) - fi)(^2 ' ■ (B.43) 

Diagrams M: 

Bl = 16ciM^ -203(^5 + 4) + 4(4 -c?i- 0^2) [(o?i-o?5)o?3' + (o?2-4)o?r']} 
+47m2[(4 - 3M^ + h/2 + i2)(M2 - 4 - hl2) 

+(de - ?>Ml + ^2/2 + h){Ml -d2- ^2/2)] , (B.44) 

Bl = 203(^2-^1) + 4(4-4 -4){[(4-4)c?r'-(^^i-4)c^3'} 

+47m2[(4 - ^Ml){s, - S2) + {h - t2)(s/2 - Ml) 

Hti/2 + ^2) {ti/2 + si)) - (^2/2 + h) {t2/2 + S2)] , (B.45) 

Bl = SciMXd-i + dji) + 2c3m[(i2 - 2Ml)dl^ + (ii - 2Ml)d^^\ 

+c'Jm[s - - 5M^ + S1 + S2 + 5/4(ii + t2)] + cl/m[di + 4 + (^i + t2)/2 
-2Ml - (si + h/2 - - M2)2rf3 1 - (^2 + ^2/2 - - M^)^^^!] , (B.46) 

Bl = 8ciM2m(4i - 4^) + 2c3m[(t2 - 2Ml)dl^ - {ti - 2Ml)d^^] 
+4M[s2 - Si + {t2 - ti)/4] + C2/m[s2 - Si + {t2 - ti)/2 
+ (si + ti/2 -m" - Mlfd^^ - (S2 + ^2/2 - - Mlfd^^] , (B.47) 
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Bl = Bl + 4(^1 - S2)/m , (B.48) 
Bl = Bl + 4(24 -di- d2) , (B.49) 
Bl = c',[{2{d, - Ml) + ^1)^3-1 - (2(^2 - Ml) + t2)dt] , (B.50) 
Bl = 4[2 - {2{d, - Ml) + ^1)^3-1 - {2{d2 - Ml) + t2)d^'] . (B.51) 

Diagram N: 

Bl^B^^ 4c4m(s2 - si)dl^ , S| = = 4c4m , Bl ^ Bf ^ -2c4(s + 3m2)d6 ^ . 

(B.52) 

Diagram O: 

S| = -B^ = 4c4m(ci3 - 4)t^5 ^ , Bl = -Bl = ic^m , Bl = -Bl = 204(4 + '^m^)d^^ . 

(B.53) 

Diagram P: 

Bl = -2Bl = 2c4{s + - 4 + 4m2[(s2 - s)^^^ + {si - s)d2^]} , (B.54) 

Bf = 2S? = 2c4{2(s2 - si) +t2-h + Am\s2 - - (si - s)c?2 ']} , (B.55) 

Bl = -2Bl = 2c4m[-8 + (2(4 - Ml) + t2)d^^ + (2(4 - Ml) + ti)rf2 ']} , (B.56) 

B^ = 2Bl = 2c4m[(2(4 - Ml) + t2)(ir' - (2(4 - Ml) + ^1)^2"']} , (B.57) 

^3^ = -2Bl ^Bj, B^^ 2Bl = Bl + 16c4m , (B.58) 

Bl = -2S| = 8c4m^(d2 ^ - ^^r^) , ^4 = = -4c4[l + 2m^(dr^ + ^^2 ^)] ■ (B.59) 
Diagram Q: 

B\ = -2Bl = 2c4[4 - 4 + 4m' (4 + 4 - 4)((^3 ' + C')] , (B.60) 

= -255* = 2c4[2(s2 - si) + t2 - ti + 4m'(4 + 4 - 4)(c?^' - d^^)] , (B.61) 

= -2^2 = Acim[-A + (m' + SM^ - s - ^3 - ti/2)d^^ + (m^ + SM^ - s - - t2/2)d^^] , 

(B.62) 

= -2Bl = Acim[{m'^ + 3Ml-s-t2-ti/2)d^^-{m^+3Ml-s-ti-t2/2)d^^]} , (B.63) 
^3^ = -2S| = -S^ , Bl = -2S| = -Bl - 16c4m , (B.64) 
Bl = -2S^ = 8c4m'(dJ^ - ^3-^) , Bl = -2S| = -4c4[l + 2m'(d3 ^ + dj^)] . (B.65) 
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Figure 1: Tree graphs with insertions from the dimension one La- 
grangian. The incoming pion is marked with an arrow. Diagrams with 
the final two pions exchanged {b c) are not shown. 
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Figure 2: Tree graphs with exactly one insertion from the dimension 
two Lagrangian. The incoming pion is marked with an arrow. The 
circlecross and the box denotes an insertion proporional to Ci,C2,c^',C3 
and C4, respectively. Diagrams with the final two pions exchanged {b <-> c) 
are not shown. 
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0.20 0.25 0.30 0.35 0.40 

[GeV] 

Figure 3: Total cross section for 7r+p — > 7r+7r+n. 
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tt' p -> it' it" n 




it" p -> it" tt° n 




0.20 0.25 0.30 0.35 0.40 



[GeV] 

Figure 6: Total cross section for 7r~p 7r°7r°n. 
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n" p -> TT° it' p 




1T' p -> IT* TT~ n 




0.17 0.18 0.19 0.20 



[GeV] 



tt" p -> TT° tt" p 
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b 




0.16 0.17 0.18 0.19 0.20 



[GeV] 

Figure 8: Threshold cross section for n'p — > n'^n'n and 
7r~p — > 7r°7r~p. 



23 



+ V + 
TT p -> TT TT" p 

5 r~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




0.16 0.17 0.18 0.19 0.20 



[GeV] 



IT* p -> it' TT* n 





2.5 


-III 


, 1 1 , 


1 


1 1 1 1 1 /I 1 
















2.0 


















( 






1.5 






















■p 
o 

4-> 


1.0 








/ y — 
/ y : 


b 


0.5 












0.0 




— 1 1 


1 


1 , , , , 1 , 



0.17 0.18 0.19 0.20 



[GeV] 



TT' p -> TT° TT° n 



3. 



b 




0.16 0.17 0.18 0.19 0.20 



[GeV] 



Figure 9: Threshold cross section for Tr+p — > 7r+7r°p, vr+p 
'K'^'K'^n and 7r~p — > 7r°7r°n. 
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Figure 10: Double differential cross section at ^/s = 1.242 GeV 
in comparison to the data of for = ooi — M^^ = 

6.5, 10.95, 15.35 MeV in the upper, middle and lower panel, respec- 
tively. 
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Figure 11: Double differential cross section at ^/s = 1.262 GeV 
in comparison to the data of for = ooi — M^^ = 

10.4, 17.5, 24.6 MeV in the upper, middle and lower panel, respec- 
tively. 
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tt' p -> IT* tt' n 
1000 I — 1 — I — 1 — 1 — , — 1 — I — 1 — 1 — I — 1 — 




Figure 12: Double differential cross section for 7T~p -K^-K^n in 
comparison to the data of at ^/s = 1.301 GeV with 0.03 < Xi < 
0.35. = ui — is the kinetic energy of the emitted tt"*" in the 
CMS. 
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Figure 13: Double differential cross section for n^p Ti^-K^n in 
comparison to the data of |^ at the incoming pion energy T^^- = 
218 MeV at xi = 0.19. T.^+ = cji — M^^ is the kinetic energy of the 
emitted 7r+ in the CMS. 
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Figure 14: Angular correlation function W for 7r~j9 TX^TX~n 
in comparison to the data of [34| for the incoming pion energy 
T^- = 284 MeV at \qi\ = 127 MeV and 62 as given in the various 
panels. 
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Figure 15: Angular correlation function W for tt p ^ tt^tt n in 
comparison to the data of |Q for the incoming pion energy T^- = 
284MeV at = 60.5MeV and \qi\ = 86MeV (left and right panels, 
respectively) at fixed 02 as given in the various panels. 
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